An exact sampling scheme for Brownian motion in the presence of a magnetic field 
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Langevin equation pertinent to diffusion limited aggregation of charged particles in the presence 
of an external magnetic field is solved exactly. The solution involves correlated random variables. 
A new scheme for exactly sampling the components of the position and velocity is proposed. 



I. INTRODUCTION 

The diffusion limited aggregation(DLA) model [l[ and 
some of its variants have been used to model pattern for- 
mation in diverse contexts such as electrochemical depo- 
sition, dissolution and erosion of porous media,gelation, 
fracture, dielectric breakdown etc 0, El 0, IE Ef These 
patterns, in general, have a fractal structure and the 
study of DLA clusters have resulted in the understanding 
of various aspects of fractals such as fractal dimension, 
multifractality and the scaling laws. The insight obtained 
from the analysis of DLA clusters has led to the applica- 
tion of this model to the study of the scaling behavior of 
two-dimensional quantum gravity. Yet a complete theo- 
retical understanding of DLA has not been achieved to 
date. Recently, a few groups 0, @ have carried out ex- 
perimental investigation of the influence of an external 
magnetic field on DLA of charged particles. They ob- 
served a bending of the branches of the DLA cluster, 
which could be attributed to the Lorentz force, as well 
as a conspicuous change in the morphology of the clus- 
ters with increase in the strength of the magnetic field. 
As a prelude to a detailed study of DLA in presence of 
an external magnetic field, we solved the corresponding 
Langevin equation exactly. It is found that the solution 
involves correlated Gaussian random variables. In order 
to carry out numerical simulation of DLA, it is necessary 
to sample these random variables, at discrete instants of 
time, consistent with the correlation relations. In this 
paper, we describe an exact sampling procedure that is 
valid for arbitrary values of all the parameters that ap- 
pear in the problem. The results of the DLA simulation 
using this prescription are in qualitative agreement with 
the experimental results and will be published elsewhere 
9]. 



II. SOLUTION OF LANGEVIN EQUATION IN 
THE PRESENCE OF MAGNETIC FIELD 

Consider the Langevin equation describing the Brown- 
ian motion of a charged particle of unit mass and charge 
q, in presence of a magnetic field B along the z-direction 



-r(i) = v(t) 
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In the above equation, if (t) — (771 (t) , 772 (i) , ^73 (^) ) is a 
Gaussian white noise of zero mean and its components 
satisfy the normalization condition 
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Since the magnetic field is oriented along the z-direction, 
it affects the motion in the x and y directions only. Defin- 
ing 

X = x + iy, V = v 1 +iv 2 , T = j + iuj, 77 = ?/i + «77 2 (4) 
we get 



— X = V 
dt 

The formal solution is given by 
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Using the above formal solution, we can immediately ob- 
tain x(t), y(t) by taking the real and imaginary parts of 
X(t). vi(t) and t>2(t) can also be obtained in a similar 
fashion. The two time correlation functions of vi(t) and 



V2(t) are given by 



(vi(t)«i(f)) = (-Wl'-^cosMt-t')) 



(9) 
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Mt)«2(t')> = («i(t)«i(f)) (io) 

(vi(t)v 2 (t')) = (J^J e -y|*-f| sin(w ( t _ f')) (ii) 

The other correlation functions can also be calculated 
easily. 

The interesting point is that unlike in the case of the 
motion without the magnetic field, now the components 
v\(t) and v 2 (t'), of the velocity, are correlated. But, the 
second moment < Vi(t) 2 >, of each of the components 
of the velocity, has the same value A/j as in the case of 
zero magnetic field. This can be easily understood once 
we recognize that < Vi(t) 2 > is twice the kinetic energy 
and that does not change by application of a magnetic 
field. However, < x(t) 2 > and < y(t) 2 > asymptotically 

rt so that the dif- 
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(i.e., as t — ► oo) tend to ]pp<- = ( 7 s + ^) <■ 
fusion coefficient D — ^i+uj^) decreases with increasing 
magnetic field. 

In the numerical simulations we will be interested in get- 
ting realizations of x(t), y(t), etc. at discrete instants of 
time tj — j * t. We have developed such a scheme that 
is valid for arbitrary values of A, ui, 7 and r. 
Define 



Xj — x(tj), yj 



Xj=Xj + iVj ( 12 ) 
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We can then write 
1 
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We can perform the numerical simulations if we can re- 
liably sample tpj and 4>j ■ since rpj and <j)j depend on 77 in 
the interval [jr, (j + l)]r] alone, and since the {rji} are 
delta-correlated, we can immediately see that (tpjip'j) and 
('Pj'P'j) vanish when j 7^ f. However, tpj and <j>j are cor- 
related. Therefore it is important to evolve a procedure 
to sample these correlated random variables. 

It is convenient to work with the following real random 
variates: 



IpjR = ^(V'j) = du T]!(jT + u) 

Jo 



du 7/2 C/t + u) 



(18) 



(19) 



jB, = $t(<f>j) = I du e 7 ^ T [cos(u(t — u))rji{jT + u) + sin(aj(r — u))r]2(jT + u) 
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[In the above equations, 3?(z) and 3(z) are the real and imaginary parts of the complex number z\. Since tpj and <j)j 
are Gaussian random variates of zero mean, the knowledge of their covariance is adequate for describing them. They 
are given below: 



((pjRCpj'R.) = 

(4>jR^3'l) = 
(4>3l*Pj'R) 



2AT8 jt j>, (V'j/V'j'l) = ZArSjj', (ipjRlpj'l) = 

^(1 - e- 2 ^)S 3 , f , (^vj) = ^(1 - e- 2 ^)^, (My!) = 



2A 



(j 2 +uj 2 ) 
2A 



[7(1 



e ' cos 



(wt)) + cue 7T sin(o;r)] Sjj 



[uj(1 — e 7T cos(wr)) — -fe 7T sin(cjT)] Sjj' 



( 7 2 +W 2) 

- ((pjR^j'l) , (<t>jl1pj'l) = ((pjRlpj'R) 



(22) 
(23) 

(24) 

(25) 
(26) 



r 



Our aim is to provide a scheme for sampling ipjR, ipji, We propose to achieve this by expressing ipjR, ipji, (pjR 
4>jH and <pji that is consistent with the above equations. and <pji as sums of n identically distributed Gaussian 
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random variables (j)ji,Pj2> •■■iPjn) with zero mean and 
unit variance in the following manner: 



1> jR = V2A^J2 
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The substitution of the above results in eauationsl22l-l2!)l 

gives 10 constraints that the parameters ai, 0i, e; and v% 
(29) should satisfy. 
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{ai}, {Pi}, {ei} and {vi) may be considered as unimodu- 
lar vectors a, 0, e* and v in an n-dimensional space. Such 
a geometrical interpretation of the above relations im- 
mediately leads to the conclusion that 4 is the minimum 
value of n, for which the above relations can be satisfied, 
see the Appendix. We consider the case n=4, henceforth. 
Since there are 16 parameters and only 10 constraints, 
there will be a 6-parameter family of solutions. We can 
choose a and consistent only with the constraints 
represented by eq. |33 and then solve for e and v. 
In what follows, we show two solutions explicitly and 
also discuss the prescription to obtain any other solution. 



A. Solution 1 

A simple choice for a and would be 



Ct\ = 1, a 2 = 0, Oi 3 = 0, «4 = 0, 
01 = 0, 02 = 1, As = 0, 4 = 



(37) 
(38) 



This choice leads to 



ei = Ci, e 2 = C 2 ,e 3 = Rcos{8), e 4 = Rsin{6), (39) 



where R = y/ (1 — C\ — Cf ) and is a free parameter. 

B. Solution 2 

Another choice for a and /3 is 



OL\ = Ct 2 = (X3 = Q?4 
Pi = -02 = -03 = 04 

We then obtained <f and v to be 



1 

2 
1 

2 



(41) 
(42) 



ei = Di + D 3 cos(6), 
e 3 = D 2 - D 3 sm(9), 
vi = D 2 + D 3 sm(9), 
v 3 = -D l + D 3 cos(9), 



e 2 = D 2 + D 3 sin(9), (43) 
e 4 = Di - D 3 cos(9) (44) 
u 2 = -Di - D 3 cos(0)j(45) 
1/4 = D 2 + -D 3 sin(0), (46) 



where D\ 



(Ci + C 2 )/2, D 2 



= (Ci - C 2 )/2, 

D3 = \/l — (C 2 + C|)/2 and is a free parameter 
which was set to be 7r/4 in our simulations. 

Any set of vectors a', 0' , e* and V obtained by simul- 

f.anpnnKlv mtn.-Hncr rotation^ a cnVpn «pt of vprt.nrs 
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a, (3, e*and v will also be a valid solution of the equations 
021 - EED It remains to be seen if a particular choice of 
these vectors is better than others for numerical simula- 
tions with finite number of realizations or not. 



III. CONCLUSIONS 

The Langevin equation for diffusion of charged parti- 
cles in presence of an external magnetic field is solved 
exactly. A scheme for sampling the position and velocity 



is given which is valid for arbitrary values of the friction 
coefficient of the medium, the strength of the magnetic 
field and the time step of evolution. The components of 
the position and velocity are represented as linear com- 
binations of four independent and identically distributed 
Gaussian random variates. A six parameter family of 
solutions is obtained for the expansion coefficients. Fur- 
ther studies are required to see if a particular choice of 
the expansion coefficients is better than all the other pos- 
sibilities. 
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The expansion coefficients oi, 0,2, b\ and b 2 can be ob- 
tained by substituting eq. in eqs. 0^1 an d E01 The 
result is 

e = Cid+C 2 (3, v = -C 2 d+C x (3 (53) 
Eq. 1531 implies 

€.e = V.v = C\ + Cf (54) 

which is not equal to 1 in general. Thus the e and v 
fail to satisfy the unimodularity condition (see eq. I48fl . 
Hence we see that it is impossible to obtain the four 
vectors a, (3, e and V satisfying the constraints eqs. 
l47!5(Jl in two dimension. 



IV. APPENDIX 

Consider four unimodular vectors a, (3, e, v satisfying 
the conditions 

a.a = l, 0.(3=1, a.f3 = (47) 

e.e = 1, v.v = l, e.v = (48) 

e.a = C u e.f3 = C 2 (49) 

v.a = -C 2 , v./3 = d (50) 

A. 2-dimensional case 

Let us first consider the 2-dimensional case. Here, the 
most general pair of vectors a and (3 satisfying eq. 03 
can be easily constructed as follows 

a = (cosQi, sinOi), (3 = {—sinGi,cos6\). (51) 

where 9\ is an arbitrary parameter. Since a and f3 are 
orthonormal, they form a basis in 2-dimensions and hence 
any other vector can be expressed as linear combinations 
of these vectors. Therefore e and V can be written as 



B. 3-dimensional case 

The most general vector a can be written down as 

ol = (sin6icos<pi,sin9isin<f>i,cos9i) (55) 

where 9\ and <pi are arbitrary parameters. The most 
general vector (3 perpendicular to a is given by 

j3 = cos(()f3 + sin(()(d x fj) (56) 

with 

fj = {cos6icos<t>x, cos9xsin<t>i, —sin9\). (57) 

and £ is a free parameter. It is obvious that the vectors 
a, (3, ax (3 form a basis in three dimensional space. Then 
e can be defined as 

e = a\d + a 2 (3 + 0,3d x f3. (58) 

Using eqs. 0^1 and then 1481 we get 

e = dd+C 2 (3 + a 3 dy. (3. (59) 
with C\ + Cf + a\ =1. Similar analysis gives 

r, — _r*„™ 4- n, r j_ h„x v r iao\ 
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with C\ + C\ + 63 =1. Substituting for e*and v, we get 

e.v — —C1C2 + CiC 2 + a 3 b 3 — a 3 b 3 (61) 

which is not zero in general. Hence it is impossible to 
obtain the four vectors a, (3, e and v satisfying the con- 



straints eqs. 14 71501 in three dimensions as well. It is, 
however, possible to obtain solutions in four dimensions, 
as explicitly demonstrated in the text. 



